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It has recently been realized that zero modes with projective non-Abelian statistics, generalizing 

the notion of Majorana bound states, may exist at the interface between a superconductor and a 
ferromagnct along the edge of a fractional topological insulator (FTI). Hero wc study two dimensional 
architectures of these non-Abelian zero modes, whose interactions are generated by the charging and 
Josephson energies of the superconductors. We derive low-energy Hamiltonians for two different 
arrays of FTIs on the plane, revealing an interesting interplay between the real-space geometry of 
the system and its topological properties. On the one hand, in a geometry where the length of the 
FTI edges is independent on the system size, the array has a topologically ordered phase, giving rise 
to a qudit toric code Hamiltonian in perturbation theory. On the other hand, in a geometry where 
the length of the edges scales with system size, we find an exact duality to an Abelian lattice gauge 
theory and no topological order. 

PACS numbers: 05.30.Pr, 74.81.Fa, 03.67.Lx, ll.15.Ha 



I. INTRODUCTION 

Systems that exhibit topological quantum order^ have 
been a focus of attention in recent years. Part of the 
interest is due to the fact that they have been proposed 
as fault-tolerant quantum memories and platforms for 
qiiantum computation^ . the paradigmatic example being 
Kitaev's toric code'^. The goal is to design architectures 
effectively governed by topologically-ordered Hamiltoni- 
ans, where qubits may be stored and manipulated in a 
physically protected way. 

Majorana zero modes, realized as superconducting 
midgap excitations in either one"* or two^^^ spatial di- 
mensions, are promising building blocks for such archi- 
tectures. Two unpaired Majorana fermions at the ends 
of a onc-dimensional (ID) superconducting wire encode 
non-locally a qubit* and, when allowed to move in a non- 
strictly ID geometry, exhibit non-Abelian statistics''"^, 
allowing to implement a non-universal set of quantum 
gates through ordered exchanges of their positions. In- 
terest in Majorana fermions has increased considerably 
in recent times, since there are now several experimen- 
tally accessible systems that may host these quasiparti- 
cles (see Rcfs. 10 and 11 for a review). A notable ex- 
ample is the edge of a two-dimensional (2D) topological 
insulator^^'^^, which hosts gapless helical (i.e. counter- 
propagating) modes, in proximity to an s-wave supercon- 
ductor (SC) and a ferromagnet (FM). The competition 
between the proximity-induced SC and FM pairing along 
the edge results in the presence of a Majorana fermion 
at each domain wall^^. 

Recently, an interesting extension of this model was 
proposed in Refs. 15-18. While the edge excitations of a 
2D topological insulator are normal electrons, it is pos- 
sible to consider instead edge quasiparticles with a frac- 
tional charge e/m, where m is an odd integer. Such gap- 
less quasiparticles appear, for example, at the edge of the 
Laughlin fractional quantum Hall states, where they arc 
described by a chiral Luttinger liquid theory^^'^°. Due 



to the absence of time-reversal symmetry, these are chi- 
ral excitations. Helical e/m quasiparticles would arise 
at the interface between two v = 1/m quantum Hall 
liquids with Lande g-factors of opposite sign or, simi- 
larly, as a Kramers doublet at the edge of a fractional 
topological insulator^^ (FTI). Also in this case, the prox- 
imity effect with superconductors and ferromagnets re- 
sults in the presence of zero modes. Since the second- 
quantization operators associated with these zero modes 
inherit a fractional exchange phase (27r)/(2rn) from the 
unperturbed edge fields, we shall refer to them as 
parafcrmionic (PF) zero modes. They are projective 
non-Abelian anyons^^'^"^, with fusion rules that generalize 
those of Majorana fermions, affording extended compu- 
tational power-^"'''^^''^^'^"''^'''. 

These superconducting zero modes realize a ID model 
with 1.2m symmetry studied by Fendley^^, which is an ex- 
tension of Kitaev's Majorana chain model** and hosts PF 
zero modes localized at the edges of the system. While 
the Kitaev chain is dual to the quantum Ising chain via 
a Jordan- Wigner transformation, the chain model is 
dual to the ID chiral Potts (p-clock) model, withp = 2m. 

Indeed, PFs as collective degrees of freedom are indeed 

well-known in statistical mechanics^^. They appear nat- 
urally in the study of the 2D p-state clock models^^"^° 
and their quantum ID counterparts^*. In lattice sys- 
tems, they arise as products of order and disorder opera- 
tors defined for self-dual systems. Moreover, PFs admit 
a description in terms of a Zp invariant conformal field 
theory^^ (CFT) featuring the PFs as primary fields. PF 
zero modes in superconducting systems are however re- 
lated only to CFTs with unit central charge (sec, for ex- 
ample, Ref. 33), arising from a bosonization description 
of the FTIs edge modes. 

In light of this body of research, it is interesting to ex- 
tend the recent works on PFs to 2D networks of su- 
perconductors. For Majorana zero modes, this question 
has already been addressed in literature^^"''® . Majorana 
lattice models can be mapped into Ising models, allow- 



2 



a) 



FTI 



SC 



SC 



FM r 

1 * 1 I * - 



FTI SC 

3*1 1 * 1 I * 



SC 



3*. 



3*N 



b) 



I * 1* 1 i * [^^ 

FTI 
N* i I * 



SC 



FTI 
N* i I * 



SC 



FTI 
N* i I * 



FM 



3*^ 



SC 



FTI 



SC 



FIG. 1: The two difFerent 2D architectures considered in this 
paper. They are composed of superconductors (SC, blue) 
and ferromagnets (FM, red) deposited on top of a 2D array 
of fractional topological insulators (FTI, grey). Z2m PF zero 
modes, marked as black stars, arise at each SC/FM interface 
along the edge. We consider two possible geometries: in panel 
(a), the FTIs extend for the whole length of the system, while 
in panel (b) the FTIs have fixed size. If we enlarged on the 
horizontal direction the system in panel (a), the number of 
FTIs would stay constant and the edge length would increase, 
while the vice versa would happen in panel (b). In the main 
text, we refer to the two architectures as the stripe and tile 
model respectively. 



1. the fractional Josephson effect, mediated by the 
tunneling of e/m quasiparticle between two differ- 
ent superconductors, and 

2. the charging energy of the superconductors. 

The fractional Josephson effect arising with J.2m PF 
zero modes has already been investigated in Refs.^^'^^, 
while to our knowledge the interplay between PF zero 
modes and Coulomb energy was not considered in previ- 
ous works. 

While both the stripe and the tile architectures give 
rise to a square lattice of PF zero modes, and de- 
spite the fact that the effective Hamiltonian contains the 
same set of local interactions in both models, we will 
show that the different geometry of the FTI edges is cru- 
cial to determine their properties, which turn out to be 
quite distinct. Indeed, since different edge geometries 
yield a different set of commutation rules for the 
PF operators and different physical constraints on the 
Hilbert space, they can determine difFerent topological 
properties. 

The paper is organized as follows: in Sec. II we derive 
the effective Hamiltonian for the stripe and tile archi- 
tectures, considering both Josephson and Coulomb ener- 
gies, and explain the physical constraints and conserva- 
tion laws specific to each array. In Sec. Ill we map the 
effective Hamiltonian into two different clock models, us- 
ing two non-isomorphic sets of PF Jordan- Wigner trans- 
formations. We analyze the phase diagram of the two 
models in Sec. IV, where we show that the tile model 
realizes a qudit toric code Hamiltonian in perturbation 
theory while the stripe model is dual to a gauge theory. 
We conclude with an outlook in Sec. V. 



ing for a description of their phase diagram. They can 
exhibit topologically ordered phases and realize the toric 
code in a perturbative regime^^. The extension of these 
analyses to ^2m PFs may reveal new Hamiltonian real- 
izations of fault-tolerant stabilizer codes^^'^® for quantum 
bits with 2to states, and hence novel platforms for quan- 
tum memories generalizing the toric code'^^. However, for 
^2m PFs the extension to 2D lattices is less immediate 
than in the Majorana case, partially because of the con- 
nection to clock models, which are less well understood 
than Ising models. 

In order to fill this gap, in this work we consider two 
distinct 2D architectures of J.2m PFs, shown in Fig. 1. 
The architectures are obtained from a pattern of super- 
conductors and ferromagnets layered on top of an array 
of 2D FTIs. The only difference between the two models 
is the geometry of the underlying FTI array. In Fig. la) 
the array is formed by long stripes of FTIs extending for 
the whole length of the system, while in Fig. lb) the 
stripes are cut in smaller pieces (or tiles) of fixed dimen- 
sion. For ease of discussion, we shall refer to these two 
architectures as the stripe and tile model respectively. 

Similarly to Refs.^^"^^, the effective Hamiltonian of the 
two models is dictated by two mesoscopic phenomena: 



II. EFFECTIVE HAMILTONIAN FOR 2D 
PARAFERMIONIC ARCHITECTURES 

In the two architectures in Fig. 1, each FM/SC in- 
terface along the edge of a FTI hosts a PF zero mode. 
Hence the stripe and tile architectures generate arrays of 
interacting J.2m PFs, which are protected by the super- 
conducting and ferromagnetic gaps and thus determine 
the low-energy physics of the system. In this section we 
derive an effective Hamiltonian in terms of J.2m PF op- 
erators for a square lattice of dimensions x Ly with 
open boundaries'*^. Generalizations to other boundary 
conditions can be easily implemented. 

Each PF zero mode is described by a second- 
quantization operator a satisfying the relations 

a^" = 1 (1) 
at = a^"-! (2) 

We can associate to a and respectively the annihila- 
tion and creation of a charge e/m on the adjacent super- 
conductor, in such a way that a Cooper pair is split in 
2m quasiparticles^^. 
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FIG. 2: A square lattice of Z2m PF zero modes of dimensions = 8 and Ly = 4. To label the PFs, we follow the convention 
established in the main text: first we order the FTI edges with an index a, and then we order the PFs along each edge with an 
index j, starting from an arbitrary origin. Each PF zero mode is then denoted as a^^jy, and all commutation rules between 
operators at different sites are fixed unambiguously. 



Furthermore, these PF operators obey unconventional 
commutation rules. Denoting two different PF operators 
with generic, ordered labels ^ and v, we have 

"m". = e-^^--/™ a,a^ (3) 
«t a, = e+--^/" a^al (4) 

where e^^, = — Ci^p = ±1 is a sign that must be fixed 
by convention. As we outline in Appendix A, Eqs. (1- 
4) can be derived from the underlying helical Luttinger 
liquid theory for the FTI edges. Note that for m = 1 
the e^ij/'s do not play any role and the equations (1-4) 
reproduce all the properties of Majorana fermions. 

In the ID case, fi and v are integers denoting the posi- 
tions of the PFs on a line. All signs are fixed by assigning 
an orientation to the line, so that e^^/ = sgn(/i — v). In 
two dimensions the ordering procedure is slightly more 
complicated and proceeds in the following way. 

1. We label each FTI edge of the system with an inte- 
ger a, thus introducing an ordering of the edges. We 
also assign a counterclockwise orientation to each 
edge a. 

2. Starting from an arbitrary origin and following the 
counterclockwise orientation, we label all ferromag- 
nets along the edge with an integer A: = 1, . . . , M 
(similarly to what was done in Ref. 15). The num- 
ber M is the total number of FMs along each FTI 
edge: note that M — 4 for the tile model, while 
M = Lj; for the stripe model (see Fig. 2). 

3. We identify the SC/FM interfaces at the left and 
the right of each FM with an integer j = 2fc — 1 or 
j — 2k respectively. 

The PFs a^,a,y are thus labelled by a composite index 
fi = {a,j}, v — {a',j'} and we fix all the conventional 
signs as 

€f,^ = sgn(a - a') -t- (5aa' sgn(j - j') (5) 



In Fig. 2 we explicitly illustrate the procedure for la- 
beling all the PFs of our square array, in both the tile 
and stripe architecture. Due to the different number and 
disposition of the FTIs, the PFs in the stripe model are 
actually distinct from (non-isomorphic to) the PFs in the 
tile model. The value of e^i, may differ for pairs of PFs in 
the same site of the square lattice, and consequently the 
set of commutation relations Eqs. (3,4) is not the same 
for the two geometries. 

In order to describe physical interactions between PFs, 
it is useful to introduce the operator 

PM. = e"-'-"/'""t"M> (6) 

defined for every given pair a^,Q!i/. For m — I, P^,^ 
represents the Z2 fermionic parity associated with two 
Majorana fermions. Here we are extending this notion 
to the J.2m symmetry of the PFs, and we shall refer to 
P^i^ as charge. From its Hermitian conjugate, 

^;.=e-— /2™ata. (7) 
we see that it is a unitary operator, 

PfJ-l^Plv = PliyPfJ-f — 1 (8) 

Moreover P^™ = 1. Thus its eigenvalues must be the 
(2m)-th roots of unity, 

A„ = e™'^/™, n = 0, 1, • • • , 2m - 1. (9) 

A pair a^, a^, can be irreducibly represented on a 2m- 
dimensional Hilbert space, with a basis given by the 
states \n) such that Pp,y|n) = A„|n). The Hilbert space 
dimension of a square lattice of PF zero modes of size 
Lx X Ly is therefore (2m)^^'^''/^. 

Now that we have set the basic algebraic rules, we 
can write down the physical ingredients of the model - 
Josephson and charging energy. These will form the ba- 
sic local bonds used to write an effective 2D Hamiltonian 
for the PFs. 
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A. Josephson energy 

Thanks to the presence of zero modes, phase-coherent 
tunneling oi e/m. quasiparticlcs may take place across the 
ferromagnetic region between adjacent superconductors 
along a common edge. The resulting Josephson effect 
is characterized by an anomalous periodicity of A-Km (in 
units of the superconducting flux quantum $o = h/2e), 
essentially because the tunneling quasiparticle's charge is 
reduced by a factor 2m with respect to the charge of a 
Cooper pair^*"^''. In other words, the anomalous period 
reflects the fact that the junction can be in 2m different 
states associated to the l-2m charge of two PFs located at 
its ends. Physically, these states are distinguished by the 
fractional spin of the fcrromagnct inside the junction'^^, 
or equivalently by the number of the fractional quasipar- 
ticles trapped in it (modulo 2m). 

Using the notation introduced in Fig. 2, the 
charge of a junction situated on edge a can be written 
as 

^'{a,2fe-l},{a,2fe} = e-^'^/^m ^fc} "{a,2fc-l} (10) 

It acts as a transfer operator, destroying one e/m charge 
inside the superconductor on one side of the junction and 
creating it on the other side. Such tunneling processes 
can be modeled by an effective Hamiltonian of the form 

ifj = (e'^'-^y^-^ "|a,2fc}«{a,2.-l} + H.C.) (11) 

Here J is the tunneling strength and 5 is the phase dif- 
ference between the two superconductors. The tunneling 
Hamiltonian splits the states of the junction in 2m energy 
branches given by 

£j,„ = -Jcos(^ + '^) (12) 
\2m m J 

with n = 0, . . . , 2m — 1. 

B. Charging energy 

Let us now consider a single superconducting island of 

our array, and let us denote with (j) and N = ~2i ^ the 
phase and number operators of this island. The presence 
of the PF zero modes becomes manifest through non- 
trivial (twisted) boundary conditions in the condensate 
ground state wave-function '^/{(f))^^, 

*(?i -h 27r) = e"« *(^) . (13) 

Here q represents the charge in units of e inside the su- 
perconductor (modulo 2e) . The spectrum of the number 
operator depends on these twisted boundary conditions, 
since its twisted eigenfunctions Xs(<?^) = e*('*+'^/^^'^/-\/27r 
satisfy 

Nxs = {2s + q)xs, seZ. (14) 



In a conventional superconductor, g = 0, we have 
periodic boundary conditions, and N counts Cooper 
pairs. In the presence of Majoranas q may assume ei- 
ther value {0, 1} giving periodic or anti-periodic bound- 
ary conditions depending on the fermion parity of the 
superconductor^^. In the presence of 1.2m PF zero modes, 
the possible values of q are extended to fractional values 

9 = {mi = i*^' m' m' • • • ' 1' • ' ' ^^^i' i^^) 

The resulting boundary conditions are twisted with pos- 
sible phases e*"'^/™, and the spectrum of TV is given by 
rational numbers with denominator m. 

Ground states with different values of q are not any- 
more degenerate if the charging energy of the supercon- 
ducting island, 

ifch = Ec{N- n^f , (16) 

is taken into account. Here Ec = e^/2C, C is the self- 
capacitance of the supcrcondiictor, and nind the charge 
(in units of e) induced on the island by nearby voltage 
gates. 

For our purposes, it is useful to separate the contri- 
bution of the fractional charges to the charging energy 
from that of the Cooper pairs. We will therefore work in 
a regime which highlights the role of the former, as done 
in Ref. 42 for Majorana fermions. If all superconducting 
islands are connected to a grounded superconductor via 
a conventional Josephson junction of energy Ej » Ec, 
the superconducting phases are pinned to their classi- 
cal minima, freezing the bosonic degree of freedom as- 
sociated with Cooper pairs. The charging energy splits 
the ground state degeneracy by inducing quantum phase 
slips. In this regime, H^h can be replaced by an effective 
Hamiltonian of the form^^ 

-ffA = -A cos(7r(jf + Trnind) • (17) 

The cosine dependence on the charge in this effec- 
tive Hamiltonian is reminiscent of the Aharonov-Casher 
effect When a (Josephson) vortex encircles a super- 
conducting island, it acquires a phase proportional to 
the charge contained in it. The energy A is exponen- 
tially small in the ratio Ej/Ec- 

Let us now write explicitly the interaction (17) in terms 
of the PF operators. We denote as qa^k the fractional 
charge trapped inside the segment of an FTI edge a be- 
tween the fc-th and (fc -|- l)-th ferromagnet. As such, it 
can be expressed as 

gi^q,,, = P^3_2fe}{a,2fc+1} = e-'^ "{a,2fc+l}"{a,2fe} (18) 

In the special case k = M, Eq. (18) has to be supple- 
mented with the boundary condition 

a{B,2M+i} = e"''"'' a{a,i} (19) 

where is the total fractional charge along the edge a. 

Eq. (19) appears naturally in the bosonization descrip- 
tion of the PFs as the boundary condition of a closed 
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edge with total charge Qs surrounding no net magnetic 
flux/^ see also Appendix A. This boundary condition 
constitutes a constraint that the physical states of the 
system must fulfill, as we will discuss more extensively in 
the final part of this Section. 

In our 2D architecture, each SC gaps out either one or 
two segments of an FTI edge, depending on whether it 
lies at the boundary of the system or in the bulk. In the 



second case, the total fractional charge q contained in it is 
the sum of two charges (/a.fe; 9a', fc' and can be expressed as 



= -P{a,fc}{a,fc+l}-P{a',/o'}{a',fe' + l}) since 



two charges operators always commute if they do 
not share a PF operator. The charging energy takes the 
form 
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I 



+ h.C.) 



on the boundary 



in the bulk 



(20) 



The total charges of the FTI edges may appear in the 
Hamiltonian (20) as additional phases, due to Eq. (19). 



C. Effective Hamiltonian 

Adding together the contributions from all islands and 
junctions, we arrive to an eff'ective Hamiltonian 



(21) 



islands 



junctions 



Each PF of the array belongs to one superconducting 
island and one junction and therefore it appears twice in 
the effective Hamiltonian, once in and once in J?j. 

Note that the effective Hamiltonian is the same for the 
stripe and the tile architectures, which share the same 
lattice, the same number of parafermions and the same 
set of local interactions. Nevertheless, the presence of 
two different sets of commutation rules for the PF opera- 
tors is enough to give the two systems markedly different 
properties. 



D. Conserved quantities and charge constraints 

The two (liffercint commutation rules between PFs in 
the stripe and tile architectures are due to the fact that 
the Hilbert spaces of the two system are constrained in 
physically different ways. To sec this, notice that the 
total charge at the edge of each FTI a must be con- 
served since no term in the Hamiltonian (21) introduces 
tunneling between different fractional topological insula- 
tors. That is. 



[e*'^*, H]=0 



(22) 



for every a. Moreover, the total charge of each FTI (edge 
plus bulk) is not only conserved but also constrained to 



be an intc;gcr multiple of the electron charge e. Thus, 
if we make the simplifying assumption that there arc no 
fractional excitations trapped in the bulk of the FTIs, 
we come to the conclusion that all ga's must be integer- 
valued. This requirement restricts the possible eigen- 
values of e*''*" to ±1, corresponding to the even or odd 
fermion parity sectors. Without loss of generality, we will 
assume that each FTI has an even number of electrons, 
Q'a = 0, ±2, ±4, ... , leading to the set of conditions 



M 



P{a,2k},{a,2k+1} = 1- 



(23) 



fe=i 



This choice amounts to restricting the twisted boundary 
conditions (19) to the periodic case. 

The constraint (23) is violated if a quasiparticle or a 
quasihole is introduced in the bulk of the FTI. Due to the 
incomprcssibility of the FTI liquid, this process is related 
to the presence of an additional flux quantum $0 = h/2e 
piercing the bulk FTI*^. Thus, we can translate the con- 
servation of electric charge on the edge of the FTI into a 
conservation of the magnetic flux threaded through the 
bulk. The latter is measured by the Aharonov-Bohm 
phase of a quasiparticle performing a counter-clockwise 
loop along the edge of the FTI. Mathematically, the 
Aharonov-Bohm phase factor is given by the string prod- 
uct Sa of the tunneling operators along such loop. 



M 



{a,2fe-l},{a,2fe} 



(24) 



fc=i 



which obeys [Sa,e* 



and has eigenvalues = e"^"/"*, n = 0, 



0, [E„H] = and (Ea)^™ = 1 



, 2m — 1. One 

can derive from Eqs. (3) and (4) the commutation rule 

-iTT/m^t y; ^25) 



*{a,j} 



""{a,.}^a 



which confirms that the operator aj^ , creating a charge 
e/m on the edge of the FTI, at the same time adds —n/m 
to the Aharonov-Bohm phase. 
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Consistently with the constraint (23), we consider as 
physical only the sector of the full Hilbert space in which 



= 1 



(26) 



so that the magnetic flux in each FTI must be a multi- 
ple of 4m<i>o. Changing the fluxes that pierce each FTI, 
we choose a different set of eigenvalues CTa and select a 
different physical sector. 

The differences between the tile and stripe model origi- 
nate from the fact that, for fixed system size, the number 
of FTIs is greater in the tile than in stripe model. This in 
turns determines the number of independent constraint 
on the Hilbert space, as reflected in the extent M of the 
product deflning the string operator in Eq. (24) - we re- 
call that il/ = 4 in the tile model and M — in the 
stripe model. The result is a different dimensionality of 
the physical sector of the Hilbert space of the two models. 
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FIG. 3: The transformation defined in Eqs. (30), (31) maps 
the PF operators living on the sites of a square lattice into 
a set of clock operators a, t defined on the hnks of the lat- 
tice occupied by a superconductor (marked in figure as black 
dots). Here, as an example, we show the positions of the clock 
operators in the case of a FTI in the tile model. The mapping 
between a\ . . . , ag and cri, ri, . . . , 0-4, r4 shown in this figure 
is explicitly written down in Eq. (37). 



III. MAPPING TO 2D QUANTUM CLOCK 
MODELS 



To highlight the differences between the two models, it 
helps to use a mapping onto quantum 2D 2m-clock mod- 
els. These models are defined on a 2D lattice where each 
site r is a 2m-level quantum system and their Hamiltoni- 
ans possess discrete l-2m local symmetries. Clock Hamil- 
tonians are defined in terms of degrees of freedom Gr and 
Tj. that satisfy 



Jim 



rt 



2 m 



(27) 
(28) 



Operators on any given site have commutation rules sim- 
ilar to those of PFs, 



G '''' "^T* f^T" 



(29) 



but operators on different sites commute. If m = 1 these 
relations are satisfied by Pauli matrices cr^,(T^. In the 
mathematical literature, the algebra describing PFs is 
known as a generalized Clifford algebra. Its representa- 
tion theory has been worked out in detail in Ref. 46, 
and from it, it is possible to infer a mapping relating the 
PF operators a and the clock operators a, r. This map- 
ping can be achieved via a parafermionic Jordan- Wigner 
transformation^^''*^''*''. 

To each couple of adjacent PFs Q;{a_2fc}: Q^{a.2fc+i} on 
the same superconductor, we associate a couple of op- 
erators <7a,kT'''a,k- Thcsc Operators therefore live on (a 
subset of) the links of the square lattice defined by the 
PFs. The mapping between PFs and clock operators is 
given by 



l<(<fe 

a{3.2k+l} = e*25S Ta^fc+i (Ja^fc+l Yl "^^'^ 

KKk 



(30) 
(31) 



Here Ka are fractional Klein factors taking care of the 
commutation rules between parafermions on different 
edges''®''*^. 



f^a l^a 



sgn(a' — a)7r/ m 



^a' ^a • 



(32) 
(33) 



These commutation rules must be compared with Eqs. 
(3,4,5). Apart from fixing the commutators, the Klein 
factors do not play a role and drop out from any 
quadratic operator considered in this paper. The bound- 
ary conditions (19) and the constraints (23), (26) are 
taken into account by setting 



(^a.M+l 

TaM+l 
M 



0-a,l 
Ta.l 



n -a = 1 

k=l 



(34) 
(35) 

(36) 



Let us write down an explicit example of the transfor- 
mation for the case M — 4, relevant for the tile architec- 
ture. In this case the relations (30) and (31), dropping 
the index a and the Klein factors for clarity, read (see 
also Fig. 3) 



(37) 



We now rewrite the Hamiltonian in terms of the clock 
operators. For the Josephson energy, Eq. (11), we obtain 



as 


cri , 


ai 


= e Tsi 


TlCTi , 




= 0-2T1 , 


as 




T2<J2Tl , 




= 173X2^ , 


"5 


= e'55S 


T30'3T2Tl , 


ae 




ar 


= 


Tia4T3T2Ti 



^l,k+l'' 



h.C. 



(38) 
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a) STRIPE MODEL 



TT + T' 
• • 



T + r' 



TT + T^T^ 



d ' a + aa 



b)TILE MODEL 





• ' T- 









tT 



a + aa' 



FIG. 4: Layout of the stripe and tile architecture in terms of clock operators a,r (black dots), living on the links of the square 
lattice occupied by a superconductor. Since clock operators at different sites commute, it is not necessary to order the FTI nor 
to assign an orientation to the FTI edges. However, notice that the clock operators for the two models live on two inequivalent 
lattices. 



while the charging energy, Eq. (20), becomes 
A , 



2 
A 



Ta,fe + h.C.) 



on the boundary 



in the bulk 



H = H, 



bulk 



boundary i 



(39) 

Note that the locahty of the interactions is preserved. At 
this point, it is useful to split the array Hamiltonian of 
Eq. (21) into bulk and boundary contributions, 

(40) 

-h.C 

(41) 
(42) 



with 
and 



junctions 



islands 
e bulk 



boundary 



A 



E 

islands 
G bdr 



see also Fig. 4. In writing Eqs. (41), (42) we have, for 
simplicity, set n-md = for all islands and, in agreement 
with our choice of the physical sector, (5 = for all junc- 
tions. Then the couplings J, A become purely real and 
all equal. 

Splitting the Hamiltonian into bulk and boundary con- 
tributions is useful for studying various boundary condi- 
tions. For simplicity, in the remainder of this paper we 
will set -ff boundary = and focus on bulk properties of the 
array, assuming the system size is large enough to jus- 
tify neglecting ^boundary Also -ffboundary = corrcspouds 
(for any system size) to the exact boundary conditions 
in case every superconducting island at the boundary of 
the array is grounded (since in that case A = at the 
boundary) . 

The Hamiltonian i/buik of Eq. (41) commutes with an 
extensive set of local operators 



3' .k' 



(43) 



associated to every bulk superconducting island s for 
both the tile and stripe models. Notice however that only 
those operators i^s that commute with the constraints de- 
scribed in the previous section are actual physical sym- 
metries: if a 7^ a', the operator (43) moves one fractional 
charge from one edge to the other, thus violating the 
charge constraint. 

The difference between the tile and the stripe archi- 
tectures can now be better appreciated, as shown in Fig. 
4. The effective Hamiltonian of H^ff Eq. (21) associated 
to either architecture is defined on the square lattice in 
terms of PFs. In contrast, iJeff is defined on inequivalent 
lattices when represented in terms of clock operators as 
in Eq. (41). The tile model Hamiltonian iJtiie is obtained 
by specializing iJbuik to a a decorated square lattice, and 
the stripe model Hamiltonian iJstripe is obtained by spe- 
cializing ifbuik to a brickwall lattice instead. 



IV. TOPOLOGICAL PHASES AND ORDERS 

The quantum phase diagram of the tile and stripe mod- 
els at zero temperature is controlled by the single parame- 
ter J/ A. In the following we will call Coulomb-dominated 
the regime A 3> J and Josephson-dominated the oppo- 
site regime J » A. In this Section we study the two 
regimes for both models, with a focus on the presence 
(or absence) of topological order. We dedicate Sec. IV A 
to the tile model and Sec. IV B to the stripe model. Let 
us summarize here the main findings. 

The Josephson-dominated regime shows no topological 
features for either model. On one hand, the ground state 
of the tile model is singly-degenerate in this limit due 
to the charge constraints. Moreover, exactly at A = 
the ground wave function reduces to a product state of 
wave functions for local four-body clusters, emphasizing 
the absence of long-range entanglement. On the other 
hand, at A = the stripe model reduces to a system of 
decoupled, one-dimensional vector Potts chains in zero 
transverse field. Hence, in the thermodynamic limit, the 
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stripe model has ferromagnetic order in the Josephson- 
dominated regime J S> A. The charge constraints do 
not suffice to select a unique ground state like for the tile 
model, but rather correlates the magnetization for pairs 
of chains. 

In the opposite Coulomb-dominated regime, and 
specifically at J = 0, both the tile and stripe mode 
show dramatically increased (relative to A = 0) ground- 
manifold degeneracy. This suggests that at least one 
phase transition separates the two regimes in the ther- 
modynamic limit for both models. 

We will show that the tile model is topologically- 
ordered in the Coulomb-dominated regime, since 

1. the degeneracy of its ground manifold depends on 
the topology of the lattice, and 

2. the model has anyonic excitations completely 
equivalent to those in the qudit toric code with l-^m 
discrete symmetry'^^ ' . 

The second point is especially noteworthy, since the tile 
model is akin but not strictly equivalent to the 1-im toric 
code by Kitaev^ nor to its generalizations^^'^^. 

Unlike the tile model, that can be defined naturally on 
a surface of arbitrary genus due to the limited extension 
of its FTIs, the stripe model fits naturally only open, 
cylindrical, or periodic (toroidal) boundary conditions. 
We will see that in the Coulomb-dominated regime its 
ground state degeneracy is not protected against local 
operators. Hence we do not consider the stripe model 
topologically ordered. We will argue nevertheless that 
the A 3> J regime is characterized by a non-local order 
parameter, that we will define using a duality mapping 
the stripe model to the l-2m lattice gauge theory^°. Thus, 
even in the absence of a topological order, the Coulomb- 
dominated phase of the stripe model can be addressed 
more generically as a topological phase. 



A. A physical realization of Z2m toric code anyons: 
the tile model 

To analyze the effective Hamiltonian for the tile archi- 
tecture, it is useful to adopt a decorated square lattice 
were each FTI sits on a site r = ie\ + je2, with 
a pair of integers, see Fig. 5. In this lattice, the array 
Hamiltonian i?buik of Eq. (41) is given by 



h.c. 

(44) 



where Q(r.r'} labels the charging energy terms of the su- 
perconducting islands, now sitting on the links of the 
square lattice between two neighboring diamonds. 




FIG. 5: Panel (a): The notation adopted for the decorated 
square lattice on which the tile model can be conveniently 
rearranged. The grey diamonds, sitting on the sites r of a 
square lattice, are FTI. To each site r there correspond four 
clock operators 0-^.,^, r^-^i, arranged counterclockwise. Blue 
links are SC, red links are FM. Panel (b): The operator Br in 
Eq. (48) is the counterclowise product of four aa^ operators 
around the same plaquette. 



Let us now consider the limit J = deep in the 
Coulomb-dominated regime. The system is then in a 
limit state where tunneling between islands is forbidden. 
Each superconductor minimizes the charging energy in 
the space of physical states specified by the charge con- 
straints: 



Q{r, 



r+ei, 



Q {r.r 



= 1. 



(46) 



These conditions allow for (2m)-fold degeneracy for each 
superconducting link, corresponding to the presence of 
the local symmetries (43). However in the sector of phys- 
ical states we must impose the constraints 



n 



1 



(47) 



(45) 



derived from Eq. (23). Nearly half of the previous states 
are then projected out, leaving a ground state manifold 
of dimension {2m)'^'^'~/'^ . This number is exact asymp- 
totically in the system size, but depends slightly on the 
boundary conditions. For example, for periodic bound- 
ary conditions the exact degeneracy of the ground mani- 
fold is (2m)i+(#SC/2)^ 

The degeneracy of the Coulomb-dominated limit at 
J = is partially lifted when weak tunneling terms are 
reintroduced, that is, we allow J 7^ 0. The Coulomb- 
dominated regime J ^ A can be treated perturbatively 
by introducing an effective low-energy Hamiltonian af- 
fecting only the ground state manifold at J = 0. We 
need to keep only those operators in the perturbative ex- 
pansion which do not couple the ground state manifold 
to the excited states. This is a standard technique^°, and 
the computation is analogous to the perturbative deriva- 
tion of the Z2 toric code Hamiltonian from Kitaev's hon- 
eycomb model^^, so we will only streamline the essential 
points. 

At first order, the perturbation cr^ ^ar,i+l creates two 
charged ±elm excitations on adjacent superconductors. 
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increasing the energy of the system by an amount G = 
2A(1 — cos7r/m). Similarly, at all odd orders we obtain 
terms that we neglect ass they do not leave the J = 
ground state manifold invariant. At second order, we ob- 
tain only terms describing the tunnelling back and forth 
of a fractional charge e/m across a single Josephson link. 
These terms renormalize the ground energy level, that is, 
they provide an energy offset to the full Hamiltonian. At 
fourth order we obtain the first relevant contribution. It 
is a plaquette operator of the form (see Fig. 5) 



Br = {<yr,i O't : 



X (7. 



'^r+ei,lO'r+ei,4 



(48) 



+62,2 y ' 



describing the tunnelling of an e/m excitation along a 
loop of four FTI edges and four superconducting islands. 
The resulting perturbative Hamiltonian reads 



frtile _ 

^pert ^ 



A 



Q{r.,r'] 



5J4 



{r.r') 



4G3 ^ 



(49) 



where we note that, in the case m = 1, the coefRcient of 
Br matches the one obtained in a similar perturbative 
expansion in Ref. 35, where equivalent plaquette opera- 
tors are obtained. 

Since the operators Br commute with the charge con- 
straints, the space of physical states for the perturbative 
Hamiltonian -ffpert is left untouched. The bond operators 
Q(ry) s-nd Br^ together with their Hermitian conjugates, 
constitute a completely commuting set of stabilizers for a 
qudit surface code"^^. This surface code protects against 
every local error that excites a ground state of iJpert into 
a state of higher energy^^. In particular, the Q^r.r') op- 
erators enforce the absence of charge excitations in the 
superconducting islands, while the Br operators enforce 
the absence of flux excitations. 

Let us note that the Hamiltonian iJpert is not exactly 
equivalent to the l-2m toric code originally discussed by 
Kitaev in Ref. 3, since the stabilizers Qi^r.r') E'-nd Br are 
not projectors. The construction of these operators is 
instead more closely related to the qudit surface codes 
introduced in Ref. 39 - although, strictly speaking, i?pert 
is not equivalent to those models as well, since it is not 
possible to canonically associate our stabilizers Q^r.r') 
and Br to vertices and faces of a two-dimensional simpli- 
cial complex. Despite these minor differences, however, 
the topological properties of these models are the same. 

The conditions 



Q {r.r 



1, 



Br 



1, 



n 



1, 



(50) 



specify the ground manifold of the system. Its degen- 
eracy can be determined from symmetry considerations. 
On the torus (periodic boundary conditions) there are 
four types of loop symmetries Ht, Ha-, Vr,Va-, defined in 
Fig. 6, associated to non-contractible loops and com- 
patible (commuting) with the charge constraints. It is 




FIG. 6: The four non trivial loop operators that define the 
ground state manifold of the tile model on a torus. The oper- 
ators Ht , Vt are defined as the product of all r along the path 
described by the two blue lines, in the order estabhshed by 
the arrows. Similarly, the operators Ha, Va are defined as the 
product of all ctct^ operators along the path given by the red 
lines. Loop operators corresponding to different paths only 
differ by a product of stabihzer operators Q{r,r') or Br- 



interesting to notice (for comparison with other sur- 
face codes in the literature) that the loop symmetries 
Hj-jVr, {Hcr,V„) are disjoint, that is, they do not have 
any clock degrees of freedom in common. As usual, any 
two loop symmetries of a given type, r or cr, associated to 
equivalent but different non-contractible loops differ only 
by a product of the stabilizer operators in (50) (or their 
hermitian conjugate). Hence, in the ground manifold, all 
these loop symmetries collapse into just four inequivalent 
ones. These form two non-commuting pairs. 



(51) 
(52) 



while [Hr,Ha] = [Vr,Va] = [Vr,Va] = [Hr,Ha] = 0. 

Since V^™ = 1 = 7?^™, it follows that each pair identi- 
hes 2m different ground states, yielding a ground state 
degeneracy of (2m)^. This is the dimension of the code 
space defined by ilpgrf 

Our stabilizer code can be adapted to a planar geome- 
try with open boundary conditions along the lines set in 
Ref. 37. In this case there will be only two non-trivial 
string operators and thus 2m ground states. In planar 
geometries with g holes, the ground-state manifold de- 
generacy increases to (2m) f. 

The Hamiltonian iJpgrt has two different types of exci- 
tations illustrated in Fig. 7: 

1. Two localized charge excitations ±e/m can be ere- 
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FIG. 7: Vortex and charge excitations in the tile model, ap- 
pearing at the ends of open string of r and ucr^ operators 
respectively. The vortices live on the plaquettes of the lat- 
tice, while the charges on the superconducting links. They 
are mutual Abelian anyons, with an exchange phase e"^''^™. 

ated on two different links by an open string of 
tunneling operator of the form S — Yii'^'^^)- The 
operator 5'^ switches the sign of the charges at the 
end of the open string. 

2. Two ±h/2e vortices are created on neighboring pla- 
quettes by one of the two operators r on the Unk 
separating the plaquettes. (The other r operator 
belonging to the same link creates the same vor- 
tices, but with opposite sign.) The vortices can be 
moved apart without further energy costs applying 
a string T — JJt of consecutive r operators sharing 
one common plaquette. 

Both charge and flux excitations are bosons when con- 
sidered separately (since different 5* operators commute 
with themselves, as well as different T operators). How- 
ever, when a charge excitation is moved in a loop 
around a flux excitation, the wavefunction will acquire 
a (Aharonov-Bohm) phase e"/™, implying that charges 
and flux excitation are mutually Abelian anyons with a 
fractional exchange phase e"/^™. This can be verified by 
computing the commutator of a pair of S and T strings 
intersecting each other. Additionally, the underlying J^2m 
symmetry allows the presence of multiple excitations of 
charge ne/m and flux nh/2e, with n = 0, . . . , 2m — 1, 
created by the n-th power of S and T operators, as in 
the usual qudit surface codes'^^. 

Let us discuss possible terms that may destroy the 
topological order. Higher orders in perturbation the- 



ory yield larger loop operators which can be decomposed 
in terms of products of Br operators and their pow- 
ers. These higher order terms commute with i?pert and 
strengthen the absence of fluxes in the plaquettes, leaving 
the ground state manifold intact. The description breaks 
down only when the perturbation order L is equal to the 
system size. At this point, the loop operators H„,V„ 
are generated in the perturbative expansion, lifting the 
ground state degeneracy by an energy 0{J^ /A^^^). 

However, we may worry about external perturbations 
of the form 

hJ2{Tr.,+Tlj, (53) 

which would break the ground state degeneracy. This 
perturbation corresponds to an external magnetic field 
in the vector Potts description and may drive a tran- 
sition from the topologically ordered phase to a topo- 
logically trivial one constituted by a condensate of the 
vortex excitations. For the Z3 toric code, this transi- 
tion was observed numerically in Ref. 52. A general 
analysis^'^'^^ of the phase diagram of Zp {p ~ 2, 3, . . .) 
Wen-Levin models^^ suggests that the transition, in the 
2+ ID transverse-field Potts universality class, is of the 
first order for any p > 3 (m > 1), thus easily detectable 
due to the discontinuity in the energy density. Our model 
however is of the vector Potts (rather than simple Potts) 
type and further investigations are required to assert the 
equivalence of the two cases for generic m. 

Finally, let us briefly discuss the Josephson-dominated 
regime. The topological order in the Coulomb-dominated 
regime of the tile model disappears when the tunneling 
terms become comparable to the charging energy. In the 
opposite extreme limit, A = 0, the FTIs decouple and the 
Hamiltonian is just the sum of the Josephson interactions 
along each diamond of the lattice in Fig 5. In particular, 
to minimize the energy, the four clock operators a for 
each FTI must be aligned and, considering the charge 
constraint (36), one obtains that the ground state of each 
FTI is constituted by an equal superposition of all the 
polarizations: 

2m-l 

J_ V ^ - kTT 

|GS)r — ^ > \crr,l = Crr,2 = O-r.3 = CTrA = e^"")- 

(54) 

Thus the total ground state is simply the product of the 
states \GS)r of all the FTIs. Due to the charge con- 
straints, it is unique independently of the topology of 
the system. 

For fixed system size, the ground state will remain non- 
degenerate also when we consider a small charging en- 
ergy contribution, A ^ J, in the Josephson-dominated 
regime. In particular the effect of applying all the charg- 
ing operators Q{r,r') in ^ closed area S is to orient all the 
clock operators a inside S. The result is the formation a 
domain wall constituted by all the links along the edge 
of S, where the clock operators are not aligned anymore. 
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The energy cost of the domain wall is proportional to 
JdS^ where dS is the number of broken Josephson links 
along the perimeter of S. Therefore, the Hamiltonian 
in the Josephson dominated regime can be seen as the 
confined phase of a loop model^^, where the loops are 
the edges of domains with different spin alignation: J 
provides a tension to the loops whereas A constitutes 
their kinetic energy. Between the topologically ordered 
Coulomb-dominated regime and the topologically trivial 
Josephson-dominated regime other phases may appear 
and the full phase diagram of the tile model deserves 
further investigations. 



B. The stripe model and the J.2m Gauge theory 

As anticipated at the end of the previous Section, the 
stripe model is naturally supported on a brick-wall lat- 
tice. It is convenient to place the clock degrees of freedom 



and 



on the sites 



{{hj) M = 0, . . . - 1, i = 0, . . . ,Lj^ - 1} 

of an Lx X Ly square lattice and distinguish between the 
two sublattices defined by the conditions (i -I- j) — even 
and {i+j) — odd. In the bulk of this geometry the stripe 
model becomes the generalization of the XXZ honeycomb 
compass model"^^ with J-2m symmetry: 



stripe 



i+j— even 
J Ly-lL^-2 

j=0 1=0 



with 



Q{i,j 



+ h.c. (55) 



(56) 



Depending on the chosen boundary condition, -ffstripe 
must be supplemented with an additional boundary term 
that we will disregard for the sake of simplicity. The 
stripe unitary operators 



SjSj+i 



n Qi 

i=0 



J=0,2,...,i, 



(57) 



represent the physical constraint on the electric charge 
of the FTIs, Eq. (36); therefore the physical states j^*) 
must satisfy: 



S,S,+,\^) = \^), j=0,2,. 



(58) 



The next task is to specify the physical symmetries of 
the stripe model. The set of non-trivial unitary operators 
that commute with -ffstripe is generated by 



1=0 



J =0,1,... 
i+ j — even. 



L,, 



1,(59) 
(60) 
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FIG. 8: Notation adopted to study the stripe model. The 
black and white sites identify clock operators a, t, distin- 
guished by a sublattice degree of freedom. Grey stripes are 
FTI, blue links are superconducors and red links ferromag- 
nets. Notice that half of the vertical links are missing, thus 
the stripe model is effectively defined on a brick-wall lattice. 
The flux threaded through a single plaquette of the lattice is 
measured by the operator B in the figure (see Eq. (62)). As in 
the tile model, flux excitations (the oriented blue circles) can 
be created by an open string of r operators: however, the ge- 
ometry constrains their movement in the horizontal direction. 



We need to specify those operators in this set that also 
commute with the charge constraints of Eq. (58). The 
symmetries Sj trivially satisfy this condition, but they 
are not all independent in the sector of physical states. 
We can keep 



J = 0, 2, • • • , Lj, - 2, 



(61) 



as an independent set of one-dimensional symmetries for 
the stripe model. As explained in the introduction to 
this section, these symmetries are spontaneously broken 
in the Josephson-dominated regime (at zero tempera- 
ture). The effective dimensional reduction displayed by 
the stripe model in this regime is intimately connected 
to the one-dimensional symmetries of Eq. (61)^^. 

The local physical symmetries of the stripe model are 
given by the minimal combination of the operators ^{i.j) 
that commute with the charge constraints and they as- 
sume the form 



t 

(i+2j) 



i+j 



even. 



(62) 



These local symmetries have an immediate interpreta- 
tion: they describe the Aharonov-Bohm phase associ- 
ated to the magnetic fluxes threading the plaquettes of 
the brick- wall lattice (see Fig. 8). 

The stripe model has no global symmetries indepen- 
dent of the lower-dimensional symmetries already dis- 
cussed. The global symmetry S1S3 . . . SLy-i, corre- 
sponding to the total charge of the system, is triv- 
ially spontaneously broken in the Josephson-dominated 
regime by the spontaneous breakdown of its one- 
dimensional constituents. Other global symmetries ap- 
pear as products of the local symmetries discussed in 
the previous paragraph, and so cannot be spontaneously 
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broken either by Elitzur's theorem^®. This suggests 
that any ordered phase of the stripe model (outside the 
the Josephson-dominated hmit) must be characterized 
in terms of a generalized, non-local order parameter^-'^. 
However this is not enough to assert that the system 
shows topological order according to our previous defi- 
nition based on the topological ground state degeneracy 
and the presence of anyonic excitation. Rather, in the ab- 
sence of a Landau local order parameter, we speak more 
generically of topological phases. 

To the purpose of comparing the tile and stripe models 
in the Coulomb-dominated regime it is useful to perform 
a perturbative analysis also of the stripe model in the 
limit A ^ J. Just as for the tile model, the first non- 
trivial term appears at the fourth order in perturbation 
theory, where the perturbative Hamiltonian of the stripe 
model becomes: 




FIG. 9: The dual lattice on which the Z^rn lattice gauge the- 
ory of Eq. (67) is defined. Clock operators now live on those 
links of the original square lattice which are marked by a blue 
dot. On this new lattice we find that in perturbation theory 
the physical interactions are given by the star and plaquette 
operators Bd and Ai defined in Eqs. (66), (69). 



^^stripe 



z+j— even 



+ h.c. (63) generators 



As in the case of the tile model, also for this Hamiltonian 
it is possible to define localized fractional charge excita- 
tions and vortex excitations as end of open strings of r 
and UG^ operators. However, for this architecture the 
vortex excitations can propagate only in the horizontal 
direction, as can be realized noticing that superconduc- 
tors of different rows share no common plaquette. 

Indeed, in the Hamiltonian -ff^'^j^ipe G^-^h row of super- 
conducting islands is decoupled from the others, and is 
characterized by a ground state degeneracy of 2m. How- 
ever, the rows between different FTIs present a non- 
physical symmetry ^ which docs not commute with the 
constraint (58). Thus, accounting for the charge con- 
straints, the overall degeneracy of the ground states in 
the physical sector is [2my^^'^^^ . Crucially, this degener- 
acy is not protected against the local symmetries ^(ij). 
Since these local operators may cause transitions between 
different ground states, the stripe model does not posses 
a proper topological order. Despite this fact, the model 
is characterized by a non-local order parameter, as we 
will discuss in the following. 

To this purpose, and more in general to investigate the 
phase diagram, it is useful to exploit the bond-algebraic 
theory of dualities^^'^^ which allows us to study the bulk 
properties of the constrained stripe Hamiltonian for large 
system size. Our strategy will be to find a duality (con- 
sistent with the constraints), mapping ifstripe to a known 
model. As shown in Refs. 31,59, quantum dualities can 
be obtained as isomorphisms of bond algebras of inter- 
actions preserving locality. In principle, we could study 
the minimal bond algebra of interactions generated by 
the bonds Q(i^^) {i+j = even) and cr(ij)CT^.^^_^.^ in iJstripe- 
However, a duality derived from this bond algebra, that 
is, an alternative local representation of these interac- 
tions, may not preserve the charge constraints of Eq. 
(58), because these constraints are not contained in this 
minimal bond algebra. Hence we consider a larger set of 



Q{i,j) J 



i = 0, • • • - 1 ; j = 0, • • • ,L„ - 2; 



' i = 0, • • • - 2 ; j = 0, • • • - 1 

(and Hermitian conjugates) for the stripe model's bond 
algebra ^stripe- That is, we are including also the bonds 
Q(ij) {i + j — odd), which are absent from the Hamil- 
tonian. Such extended bond algebra does contain the 
charge constraints in Eq. (58); hence a duality for .4stripe 
maps these constraints in a well defined fashion either 
to the identity operator (in which case the duality solves 
the constraints'^^) or to constraints of the dual model. 

The characterization of .Astripe in terms of relations 
among its generators reveals the following dual repre- 
sentation of the algebra of interactions. 



Q 



(i.j) 



B. 



(64) 
(65) 



with 



if i = 0, 

^ otherwise. 

(66) 

Then the dual Hamiltonian, Hq = 'I'd(-ff stripe), reads 



rA jLy-lL^-2 



■i+j— even 



j=0 i=0 



h.c. 

(67) 

and it is unitarily equivalent '^-'^ to the stripe model. 

Up to boundary terms, i.e. the incomplete plaquettes 
B{o,j) (j — 0,2, . . . ,Ly — 2), and a redefinition of the 
lattice that places the clock degrees of freedom on links 
rather than sites, we recognize Hq as the Hamiltonian 
of the Z2„i lattice gauge theory studied in connection to 
the problem of confinement in QCD^" (see Fig. 9). The 
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local symmetries of the stripe model, Eq. (62), 

map under duality to 



A. 



with 



A. 



(68) 



(69) 



As is guaranteed by the formalism, the unitary oper- 
ators Aci(j j), (i + j = even), commute with the dual 
Hamiltonianffc • They correspond to the gauge symme- 
tries of the J-2m gauge theory and they have the inter- 
pretation of measuring the local density of external J-2m 
charge. Hence our duality maps the magnetic fluxes in 
the stripe model, as described by the Aharonov-Bohm 
operators Sj-j j) in Eq. (62), to external electric 
charges in the gauge theory. 

At this point we can exploit Eq. (68) to compute the 
dual representation H^^^^^ = $d(-ffstripe) °f the pertur- 
bative Hamiltonian of Eq. (63), 



rrpert 



E 

I + J— even 



5J4 
4G3 



Ad 



-l-h.c. . (70) 



Remarkably, this is the Hamiltonian for the qudit toric 
code model^^'^^. To clarify the notation, notice that due 
to our definition of the plaquctte operator Eq. (66), the 
two operators A^ j) , f^i j) share one vertical link of the 
lattice, with A^ (^i jj to the right and B^ to the left of 
that link (see Fig. 9). This duality, however, is a non- 
local transformation with respect to the clock operators a 
and r. Thus, even if the spectrum of H'^^J^^^ is equivalent 
to the Z2m toric code, the stripe model in the Coulomb- 
dominated regime does not present topological order. 

To assert that the phase diagram of the gauge theory 
and the stripe model are connected by the duality ^d, 
we need to investigate the effect of the duality mapping 
on the charge constraints of Eq. (58). Remarkably, the 
charge constraints are holographic^^, since they map to 
boundary constraints for Hq, 



SjS- 



'j+i 



'(L,-lj)'"(L,-lj + l)' 



(71) 



0,2,. 



2. Then the physical states 



for 

|^)d = ^'dl^') for Hq, seen as a dual representation of 
the stripe model, are characterized by the condition 



a 



(ix-lj)'^(Lx-lj + l) 



*) 



l*>d 



(72) 



for j = 0, 2, 



, Ly — 2, and not by the standard con- 



dition of gauge invariance [that is, invariance under the 
A{i.j), {i + j — even)].^" Despite this difference we will 
argue in the following that the stripe model and the Z2m 
lattice gauge theory share the same phase diagram. 

The dual gauge theory Hq presents unusual open 
boundary conditions. Since the charge constraints are 
holographic, this is in perspective required to guarantee 




FIG. 10: The generalized order parameter (string tension) in 
both the dual and original lattice. 



that the dual charge constraints supported on the bound- 
ary commute with the dual Hamiltonian Hq (just as the 
charge constraints commute with -ffstripe)- 

However, the standard view that boundary conditions 
do not affect the phase diagram in the thermodynamic 
limit suggests in this case that the dual charge constraints 
do not affect the phase diagram of Hq, which must then 
coincide with the standard phase diagram of the 
gauge theory. This view is strengthened by the fact that 
the dual charge constraints commute with the gauge sym- 
metries of Hq , and so the ground state of Hq will belong 
to the gauge-invariant sector even after the condition Eq. 
(72) is imposed. Finally this implies that the stripe model 
presents the same phase diagram independently on the 
choice of the charge of each FTI edge in the constraint 
(23), thus in all the different physical sectors. 

In view of these considerations, we can argue that the 
stripe model shares the phase diagram of the 1.2m gauge 
theory as described in Ref. 60 (and references therein). 
It follows that there is indeed one second order phase 
transition separating the Coulomb-dominated from the 
Josephson-dominated regime. In the gauge-theory lan- 
guage this transition is understood as a confinement- 
deconfinement transition. In particular, the Coulomb- 
dominated regime of the stripe model is dual to to the de- 
confined phase of the gauge theory, while the Josephson- 
dominated regime corresponds to the confined phase. 

The phases of a gauge theory cannot be distinguished 
by a Landau order parameter^*. However, for the 
gauge theory dual to the stripe model, there exists a gen- 
eralized order parameter, the so called string tensionP^, 
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that is non-zero in the confining phase and vanishes con- 
tinuously, but non-analitically at the transition point. 
The string tension is the expectation value of a string 
of r's in the tim gauge theory, which correspond to an 
open string of tunneling operators cra^ in the stripe model 
(see Fig. 10, analogously to the string operators creating 
charge excitations in the tile model. The ground state 
expectation value of such string falls continuously, but 
non-analytically to zero at the second order phase transi- 
tion separating the Josephson-dominated regime (where 
it is different from zero) from the Coulomb-dominated 
regime. On the gauge theory side of the duality, the two 
phases can be distinguished also by a different scaling 
of the expectation value of the Wilson loops, which map 
to sets of r operators in compact regions in the stripe 
model. 



V. CONCLUSIONS AND OUTLOOK 

In summary, we have studied two-dimensional arrays of 
interacting parafermionic zero modes. Such exotic states 
form along the edge of fractional topological insulators, 
at the domain walls between proximity-induced super- 
conducting and ferromagnetic pairing. The dynamic of 
these zero modes is dictated by two competing effects, the 
charging energy of each superconducting island and the 
fractional Josephson tunneling of quasiparticles between 
different islands. 

The underlying fractional edge modes, which are orig- 
inally described by a helical Luttinger liquid theory, de- 
termine crucially the possible lattice geometries and the 
physical constraints of these parafermionic systems. We 
have analyzed two possible architectures, the tile and the 
stripe model. They differ mainly for the fact that in the 
former the length of the edges is constant, while in the 
latter it scales with the total size of the architecture. 
We have discussed how this feature gives rise to different 
physics, despite the fact that the models are character- 
ized by the same lattice of parafermions and the same 
local interactions. 

The difference is appreciated by exploiting a Jordan- 
Wigner transformation mapping the parafermionic oper- 
ators into clock operators. Through this transformation 
the tile model is described by a Hamiltonian on a dec- 
orated square lattice whereas the stripe model becomes 
a compass model, with l-im symmetry, on a brick-wall 



lattice. 

The tile model presents, at least in perturbation the- 
ory, the same topological order of the surface codes char- 
acterized by a l-im symmetry: it shows the same topolog- 
ical degeneracy of the ground state and the same anyonic 
excitations. Thus the system we described is a possible 
physical candidate to the realization of qudit surface code 
Hamiltonians. It is known that the ground state degen- 
eracy of these systems suffers from thermal fragility®^ 
However, we note that the intrinsic noise due to the pres- 
ence of induced charges on the superconducting islands 
could help localize the anyonic excitations of the system, 
and thus to protect the information which may be en- 
coded in the ground states. ^''^^^ 

The stripe model provides instead a physical realiza- 
tion of the Tim lattice gauge theory, a toy model often 
exploited to study confinement-related problems in lat- 
tice field theory. The duality mapping between the stripe 
model and the lattice gaiigc theory is non-local in terms 
of the single clock degrees of freedom, but it is local in 
terms of the interactions. Unlike in the tile model, a 
toric code Hamiltonian can only be retrieved in the dual 
theory, where the operators are non-local. Interestingly 
enough the physical charge constraints of the FTl edges 
maps to an holographic constraint in the gauge theory 
which affects only boundary terms. 

To conclude, our work addressed the problem of find- 
ing topologically-ordered phases in the phase diagram of 
these two-dimensional collections of topological defects. 
Comparing the results obtained for the two architectures, 
we can see that that it is not only the nature of the inter- 
actions between the defects (in this case, the PF zero 
modes) that determines the presence of topological order, 
but also the intrinsic geometry of the topological phases 
originally generating the defects (in this case, the edges 
of the fractional topological insulators). Understanding 
the interplay between this two aspects is crucial to design 
topologically-ordered architectures. 
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Appendix A: Description of the system through 
bosonization 



In this Appendix wc summarize the main features of 
the bosonization description of our system and we pro- 
vide an expression in terms of massless bosonic fields of 
the parafermion operators a and thus of the related in- 
teraction terms. We follow the approach in Refs. 15 and 
16, where more details can be found. 

In absence of the interactions provided by the super- 
conducting islands and the ferromagnetic insulators, the 
edge of the FTIs defining our systems, or, equivalently, 
the double edges of juxtaposed fractional quantum Hall 
layers with opposite polarization, can be described in 
terms of the Luttinger liquid Hamitonian proposed by 
Weni^'^°. In particular the massless edge modes are de- 
scribed by the following Hamiltonian: 



Hn = 



27r 



da; 



(Al) 



where v is the speed of the two counterpropagating modes 
and (fi, are dual massless bosonic fields obeying the com- 
mutation relation: 



m 



{X2 - xi) (A2) 



where O is the Heaviside step function. In particular for 
each FTI edge a it is possible to define two chiral bosonic 
fields 



'Pl/r {x±t,a) = (p {x, t,a)TO {x, t, a) 



(A3) 
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in such a way that the left and right fermionic modes, 
with opposite spin polarization, are defined by the op- 
erators ^i^/ji{x,t,a) = o""'''-f-/"(^*''''^) where r/s are 
fermionic Klein factors. The charge density associated 
with each edge is p = dxO/w, thus, in a closed edge with 
length jC, the total charge of the edge is related to the 
boundary conditions of the 6 field: 

TTqtot{a)=0{x + £.,a)-e{x,a) (A4) 

and analogous conditions relate the field (p with the spin 
density 

For each edge the interaction terms corresponding to 
the proximity induced superconducting coupling and the 
backscattering give rise to the interaction Hamiltonian: 

Hi oc J dx[— gs{x) cos {2m(fi) — gF{x) cos {2m6)] (A5) 

where gs and gp describe respectively the position depen- 
dence of the induced superconducting and ferromagnetic 
couplings. 

By selecting a position in the bulk of either a super- 
conducting or a ferromagnetic segment of the edge, if 
the couplings g are strong enough, one can consider re- 
spectively the fields ip and 9 as pinned to the semiclas- 
sical minima ipf^, 9^. = 0, n/m,, 27r/m, . . . , (2m — l)7r/m. 
Adopting this approximation and considering the limit 
of sharp transitions between superconducting and ferro- 
magnetic regions, the parafermion operators can be writ- 
ten as: 

a2fc-i,a = K^e'^"^"''-'"-^ (A6) 
a2fe,a = K^e'^'fi>'+'■'^-'^''^'^^ (A7) 

where fc = 1, . . . , M labels the ferromagnets and the su- 
perconductors along the edge a and the tile and stripe 
models are characterized respectively by M = 4 and 



M = Lx- The fractional Klein factors Kg enforce the 
correct commutation rules (3) and (4) . This definition of 
the parafermionic modes is not unique (see Refs. 15 and 
16 for more detail) but it suSices to our purposes. Finally, 
for a complete description of the system, it is necessary 
to take into account the correct boundary c;onditions. 

Through this definition of the parafermionic operators 
it is easy to verify that the tunneling operators assume 
the form c-iivk+i.^-Vk,.) = P{a,2fc-i}.{a.2fc} and thus wo 
recover the usual form for the fractional Josephson inter- 
action (11): 

Hj = - J cos ^<^fe+l,a - V'fe.a - (A8) 

Moreover the tunneling string operator Eg 
defined in (24) for the two models become 
Sa = cxp (— i ((/^M+i.a ^ </'i.a)), cuphasizing the re- 
lation between the boundary conditions of the field ip 
and the magnetic flux enclosed by the FTI edges. The 
boundary condition (19) assumes a natural form in the 
bosonized description due to the boundary relation (A4) 
since: 

«{a,2M+i} = K3e'('^"+i--^-+i-) = e--«» a{a,l}Sl 

(A9) 

Once we apply the parafermionic Jordan- Wigner trans- 
formation (30,31) to map the system in a quantum clock 
model, the previous boundary conditions are translated 
in the following relations: 

M 

fe=i 

= e-iCVM+l.a-V'l.a) = <jt^^^CTl (All) 

which generalize the boundary conditions (34,35,36). 



